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NOMENCLATURE 


A Unreefed swept area of rotor 

A' Swept area of rotor when it is reefed 

A^ Stream tube area far away from tiie disc 

A Stream tube area at the wake 

W 

a Axial interference parameter at the disc 

b Axial interference parameter at the wake 

( Quarter blade length ) 

C Chord length 

Section drag coefficient 
Section lift coefficient 

Normal force coefficient acting normal to the chord 

Tangential force coefficient acting along the chord 

Cp Power coefficient of wind turbine 

Cp' Power coefficient of wind turbine based on V 

Thrust coefficient of wind turbine 

C^' Thrxist coefficient of wind turbine based on V 

Cq Torque coefficient 

dD Elemental drag 

dL Elemental lift 

dQ Elemental torque 

dT Elemental thrust 

E Generator efficiency 

E. Drive train efficiency 

trans 

G Ratio of quarter blade length to strut length 

k Theoretical = 16/27 

k^ = ( A/R^) 

L Strut length 



Number of blades 


Free stream static pressure 
Pressure acting in front of the disc 
Pressure acting on the back side of the disc 
Mean power produced by the turbine 
Mean power produced by the half blade 
Mean torque on the rotor 
Unreefed radius of the turbine 
Radius of the turbine when it is reefed 
Mean thrust acting on the turbine 
Relative velocity of the blade in the plane of the 
blade section 

Relative velocity of the blade in the horizontal plane 

Free stream wind speed 

Wind speed at the disc 

Angle of attack 

Arrowhead angle 

Azimuthal angle 

Reef angle 

Free stream air density 
Angular velocity of the turbine 
Strut angle 
Strut angle for F = 0° 

Solidity - NG/R 

Q R 

Tip speed ratio based on — ) 

Q R 

Tip speed ratio based on V ( = ) 

Rated angular velocity 
Half blade length 



Blade chord length 
M Mass of the half blade 

D 

cg^ C.G. Location of the half blade 
W Width of the end plate 

t Thickness of the end plate 

Length of the end plate 
M Mass of the end plate 

W width of the sleeve 

S 1 

t Thickness of the sleeve 

S 1 

h^ ^ Height of bearing sleeve or pad 
bear_s Space between the bearings 

Location of the strut hinge 

X Location of the near end of X - arm ( Unreefed ) 

X 

M Mass of the X - arm 

X 

(p ^ Azimuthal angle at t = 0 

Velocity at the disc after jump 
oc^ Angle of attack at t = 0 

ttj Angle of attack after jump in velocity at the disc 

Aerodynamic torque at t = 0 

Qj Aerodynamic torque after jump 

Q Generator torque 

g 

Q Aerodynamic torque produced by the half blade 

hb 

Radial velocity of the cross-arm 
u Angular velocity 

L Strut length 

.'S', 

F Net outward force in the cross-arm 

■' ■ C X 

Force in the upper hinge (Fig. 6. 7) 

F^' Force in the lower hinge (Fig. 6. 7) 

r C.G. location of the cross-arm (unreefed case) (Fig. 6. 6) 



r' C.G. location of the cross-arm (reefed case) (Fig. 6. 7) 

h^ Height of the upper hinge from sleeve 

strut angle 
F Strut force 
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ABSTRACT 


A Variable Geometry Vertical Axis Wind Turbine with a novel 
control system proposed by Dr. Ghosh in 1990, called as 
centrifugal reef control system, has been chosen for the present 
study for having certain advantages over other types of wind 
turbines. The aerodynamic design parameters such as radius of the 
turbine, blade length and chord have been chosen so as to generate 
approximately 10 kW of power at a wind speed of 10 m/s. Templin' s 
procedure, which avoids iteration, has been followed in the context 
of Rankine Froude Momentum theory to relate the induced velocity 
at actuator disc to free stream velocity. To solve equations 
numerically, Newton - Cotes technique has been used and a FORTRAN 
code has been developed for this purpose. For simplicity, the 
induced velocity is assumed to be uniform over the entire swept 
area of rotor. Blade element analysis has been used to derive the 
aerodynamic performance characteristics ( i.e. power, thrust and 
torque coefficients ) of the wind turbine. The aerodynamic 
characteristics have also been studied for reefed configurations. 
The effect of solidity on power coefficient has been presented. 
The power coefficient for different reefed conditions has been 
compared with other available theories and experiments. 
Aerodynamic forces acting on the blades have been estimated for 
the unreefed condition when they are at their maxima. A summary of 
the work done on VGVAWT at IIT Kanpur has been presented in the 


end . 



INTRODUCTION 


Winds result from the fact that the earth' s equatorial 
regions receive more solar energy than the polar regions, and this 
sets up large scale convection currents in the atmosphere. 
Meteorologists estimate that about 1% of the incoming solar : 
radiation is converted into wind energy. The solar energy received 
by the earth in just 10 days has an energy content nearly equal to 
the world's entire fossil fuel reserves ( coal, oil & gas ) . This 
means that the wind resource is - extremely large. As the fossil 

fuels are exhaustible, sufficient power will not be produced using 

) ; 

these conventional sources of energy in the near future. 
Consequently, we have to go in for alternative sources of energy 
which are renewable like Solar, Wind and Biomass energy etc. 

It has been realized by now that wind is a clean (non - ^ 
polluting) and renewable source of energy from which electrical 
and mechanical power may be generated. Right from the beginning of 
civilization, people used to harness the wind energy for useful 
purposes like pumping water, grinding grains, sailing ships etc. 
Over the last 2000 years or so, different types of wind machines t 
have been invented. Wind energy research and development work are ] 
being done in U.K., USA, China, Japan, India and many other ; 
countries . ^ 

In India, the Government has also realized the potential of | 
renewable resources of energy and it has taken various steps in 

this regard. In 1987, Government of India set up Indian Renewable 

' ■ ' ' '' ; ■' '■ : . ■ ,, . ' ■ ' " ■■ ' ■ ' ' ' ,, . ■ 

Energy Development Agency ( IREDA } to finance commercially viable 

projects . I 

Rapid progress is being made in the development and j 
demonstration of wind energy. Till March 1995, a capacity of about 
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275 MW of wind energy generation has been established; including 
private sector projects [26] . The Govt, has extended a package of 
incentives for wind power generation. This includes loans from 
IREDA at concessional interest rates, tax holiday for 5 years, 
exemption from excise duty and exemption or reduced rate for 
custom duty on imports of wind turbines. IREDA has initiated the 
'' India Renewable Resources Development Project '' consisting of 
three projects having ' Wind - farm Development ' as one of the 
projects. Till March 1995, IREDA had sanctioned nearly 400 
renewable energy projects for about Rs . 3 60 crores for manufacture 
and installation of new and renewable sources of energy systems 
and devices. 

Department of Non-conventional Energy Sources ( DNES ) has 
started its programmes at various places all over India like 
Tuticorin (Tamilnadu) , Mandvi and Okha (Gujarat) and Puri (Orissa) 
etc. in cooperation with various State Electricity Boards in 1986 
[12]. The Govt, of India has announced various policy measures to 
encourage direct foreign investments and collaborations. 

A Variable Geometry Vertical Axis Wind Turbine (VGVAWT) has 
been chosen for the present analysis for having its own 
distinguished advantages over Horizontal Axis Wind Turbine (HAWT) 
from design point of view though it has a few disadvantages also. 
A critical and comparative assessment of both VAWT and HAWT has 
been made and VGVAWT with centrifugal reef control (proposed by 
Dr. Ghosh) has been chosen. Wind Energy survey in India 
identifies many wind rich areas all over the country. The average 
annual mean wind speed is found to vary from 6 m/s to 10 m/s in 
places like Tuticorin, Okha etc. [12] . 



CHAPTER 


1 


Literature survey : 

In 1930 Darrieus first developed a VAWT which utilizes the 
lift force. Before his period, people used to try with lift; 
devices and drag devices separately. Darrieus used a blade with 
Troposkien shape, which can be explained as following : 

When a rope is rotated about a vertical axis at a particular 
r.p.m. , it takes a special shape called as 'Troposkien shape'. 
The centrifugal force acting at each point on such curve will be 
balanced by the tensile force of the rope. Hence, it is 
advantageous to use blades of this shape. But making Troposkien 
shape without bending a straight blade (i.e. without introducing 
pre - stresses) is very difficult. 

The power produced by the turbine depends on the wind 
speed. Most of the wind turbines are located in coastal areas: 
since they have high wind power potential throughout the year. The 
coastal areas are often prone to cyclones and hence a wind turbine 
with the straight blades with fixed geometry may experience severe 
loadings. Musgrove [1976] suggested a variable geometry vertical 
axis wind turbine and a 25 m. diameter, 160 kW prototype machine 
was installed in Carmarthen Bay. Reefing was achieved 
hydraulically through a computer controlled system. The experience 
and the data gathered on this machine was utilized in constructing 
a 17m diameter 100 kW variable geometry VAWT on the Isles of 
Scilly (Mays et al, 1987). In late 1987, there was a failure in 
the reefing control mechanism which led to an overspeed situation 
causing damage to the rotor (Mays 1989, Lindley 1988/89). Data' 
gathered from the 25m wind turbine at Carmarthen Bay has shown 
that a stall regulated fixed geometry rotor gives less severe; 
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loading than had been anticipated (Mays 1989). hater the wind 
turbine was replaced by a stall regulated fixed geometry rotor 
where beyond certain wind speed, the rotor is held at a constant 
r.p.m. by the grid causing the blades to stall and thus the wind 
turbine is controlled. 

The variable geometry concept is a novel idea to get a 
constant rated power and maintain a constant r.p.m. The angular 
speed should be maintained constant in the case of electricity 
generation to achieve a steady frequency, for ex. 50 Hz, in the 
case of electric power generation for the grid. 

On comparing VAWT & HAWT, we see that both have advantages 
as well as disadvantages. The VAWT is able to extract energy from 
the wind independent of the wind direction, i.e., VAWTs are 
omnidirectional so that one can avoid the cost and complexity of a ! 
yaw orientation sub - system. Due to the fact that the rotor has a 
fixed pitch and low solidity, VAWT is not self starting. One of 
the attractions of the vertical axis configuration is that it 
eliminates the cyclically varying stresses which the blades of 
HAWT experience, due to their weight, as they rotate. These 
cyclically varying gravity stresses become progressively more 
dominant as the overall turbine size increases, and imposes a 
theoretical limit on the size of the horizontal axis rotor. 
Although VAWT blades experience cyclically varying aerodynamic 
loads, the resulting stresses remain substantially constant as the 
turbine size is increased. 

'I 

The novel idea of reefing (i.e. deflecting the blade w.r.t. 
a vertical axis to resemble an arrowhead) was achieved (Anderson 
et. al.,1987) by a hydraulic system which senses the over speed of | 
the wind or wind turbine and allows the blade to reef. The | 



conditions and rotor was replaced by a fixed geometry stall 
regulated one. Dr. K. Ghosh (1990) suggested a novel method of 
reef control viz. centrifugal reef control technique. He argues 
that since the radial component of the aerodynamic forces acting 
on a blade is outward during one half of the rotation and 'inward 
during the other half, its contribution in assisting the 
centrifugal forces on an average can be ignored. He showed that 
with a realistic mass distribution among various members (i.e. 
blade, cross - arm etc.) , the centrifugal force varies fairly 
linearly with distance of the C.G. from the axis of rotation. This 
force, therefore, can be balanced against a linear spring [Fig.^.^ . 
To limit the to and fro motion due to aerodynamic force changing 
direction, dampers are necessary [Fig. 3.2]. He developed a method 
to find the centrifugal force and relation between the reef angle 
and the change of radius with reef angles. Hence it is possible to 
assess the stiffness of the spring which is capable of 
withstanding the centrifugal force for reefing angle from 0 to 90 
degrees from his analysis. His analysis also shows that the 
centrifugal reef control is cheap and quite reliable.. Moreover, 
the control system is simple to design. 

Design of such wind turbine involves thrust and power 
calculations based on geometry for various wind speeds and rpms. 
Templin [1974] deviced a methodology to find out the average power 
and thrust coefficients which avoids an iterative procedure . The 
above procedure was adopted by Shankar [1976] and he calculated 
the performance of a Darrieus wind turbine assuming both uniform 
induced velocity and non - uniform induced velocity approaches . 

Rankine Froude momentum theory further 
developed by Glauert [1963] , wherein the wind turbine is assumed 
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and interference factor. This theory is generally user] by any wind 
turbine designer for preliminary design calculations. But one 
should realize that this theory will not be adequate to analyze 
what exactly happens to the flow that passes through a vertical 
axis wind turbine. The actuator cylinder flow theory discussed by 
Freris [1990] , predicts the flow velocities just upstream and 
downstream of the turbine . 

In 1992, Selvaraj , an M.Tech. student in Aerospace Engg. 
Department at I.I.T. Kanpur worked on VGVAWT under the guidance of 
Dr. Ghosh . His formulation of the problem and computer program 
have been rechecked now (1996) . There were some errors in his 
computer program and those have been removed now. His graphs have 
been rechecked and some modifications in his formulation have been 
made as suggested by Dr. Ghosh. The computer program has been 
modified accordingly. 



CHAPTER 


2 


2.1 Aerodynamic Modelling : 

2.1.1 Rankine Froude Momentum Theory : 

This theory is based upon following assumptions : 

1) The wind is steady and homogeneous 

2) The flow is inviscid and incompressible 

3) The flow is one dimensional and has no rotational velocity 

4) The wind turbine is considered as a simple actuator disc 

5) The wind has constant velocity over the disc and over the 
cross-section of the ultimate wake. 

The flow past a wind turbine is as shown in Fig. 2.1 . The free 

stream wind, V is slowed down to V at the disc and it is further 
00 

decelerated to u. During this process, wind machine extracts 
energy from the wind. The fluid which has passed or will pass 
through the wind turbine disc forms a cylindrical duct whose 

cross-sectional area increases from A far in front of the disc to 

00 

A at the wind turbine disc and finally to A at the ultimate wake 

W 

(Fig .2.1) . 

To proceed further with the development of the theory, it is 
necessary , to consider the form of the momentum equation as applied 
to individual annular elements of the wind turbine , i . e. , to the 
fluid which passes through the turbine between the radii r and 
(r+dr) . 

The system of forces acting on the control volume of fluid 
bounded at the ends by the two sections, A^ and A^ far upstream 
and downstream the turbine Goinprlsen of thrust T, iJie pressure 



force P (A - A ) over the ends of the control volume and an axial 

00 W 00 

force X due to the pressure on the lateral boundary of control 
volume. Then the momentum equation for the control volume of fluid 
can now be expressed as 

T-X + P(A - A) = Xp u{V - u) dA ....(2.1) 

00 M 00 00 00 W 

where the integral extends over the cross-section of the wake and 
the velocity u may have different values for different annular 
elements. 

It has been shown in Ref. [10] that the axial pressure 

force X on the lateral boundary of the control volume of fluid is 
equal to the pressure force P (A - A ) acting at both ends. 

00 w 00 

Therefore the Eqn.(2.1) becomes 

T = X p u (V - u) dA .... (2.2) 

Due to the assumption no. (5) , the above eqn. may be written in a 

very simple form as 

T = p A u (V - u) (2.3) 

00 w 00 

By continuity, we have - 

m = pA V = pAV=p Au ...(2.4) 

*000000 00 *^00w 

Introducing Eqn. (2.4) into Eqn . (2 . 3 ) , we get - 

T = m (V - u) = p A V (V - u) .... (2.5) 

Applying Bernoulli's equation between the sections 0 and 3, and 2 


and 1 (Fig .2-1) , we have - 

P -t- (1/2) p + (1/2) p .....(2.6 a) 

00 00 00 00 

P" + (1/2) p = P + (1/2) p u^ (2.6 b) 

00 00 CO 


Also the thrust acting on the air must be equal to the thrust 
acting on the wind turbine which is nothing but the area times the 
pressure difference acting on either side of the turbine, 
i.e. , ; : 

T = A (P"^- P~) .... (2.7) 


P'^ - P = (1/2) p (V^ - u^) ...(2.8) 

Therefore Eqn. (2.7) becomes 

T = (1/2) p A (V^ - u^) ..,.(2.9) 

tX) CO 

Comparing Eqns.(2.3) and (2.8), it is seen that 

V = (V + u) /2 ... (2 .10) 

00 

Thus the velocity at the disc is the average of the upstream and 
downstream velocities. 

Defining axial interference parameters, a and b as the 
fractional decreases in wind velocity at the rotor plane and wake, 
represented by 



V = V 

(1 - 

a) 


00 

and 

u = V 

(1 - 

b) 


( 2 . 11 ) 

....( 2 . 12 ) 


00 

respectively. From Eqn. (2.10), it is obvious that b = 2a. From 
Eqns.{2.11) and (2.12), we can infer that when a = 0, no power is 
absorbed, whereas for a = 0.5 (i.e. b = 1) , the wind has zero 
velocity at the wake. When a = 1, the flow has zero velocity right 
at the disc and hence the streamlines break away from the turbine - 
Because of this flow behaviour at the disc, the stream tube theory 


i.e. the Rankine Froude's momentum theory can not be used in this 
condition of operation. 

The energy removed by the rotor per unit time is 

P = (l/2)p A V - (l/2)p u^A V 

00 00 00 


i . e . , 

P = (1/2) p A V (V + u) (V - u) ...(2.13) 

00 00 00 

Substituting for V from Eqn. (2.11) and u from Eqn. (2.12), the 
above Eqn. becomes 

P = (1/2) p A (4 a) (1 - a)^ .....(2.14) 

00 00 

Defining the power coefficient Gp as 

= P 

^ (1/2) p A 

*^00 00 


(2.15) 
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Then Cp , while substituting for P from Eqn . (2.14) becomes 

Cp = 4a (1 - a) ^ (2.16) 

Similarly, we can define the thrust coefficient as 

(1/2) p A (2.17) 

CO CXJ 

Introducing the Eqns.(2.11) and (2.12) in (2.9) to get T , and 

then from Eqn. (2.17) becomes 

= 4a (1 -a) (2.18) 

By maximizing Cp , it can be easily shown that the max. power 

can be extracted when a = 1/3. The corresponding value of Cp is, 

(C^ ) = 16/27 ....(2.19) 

The wind velocities at the max. power condition are 

V = 2/3 V (2 .20 a) 

00 

u = 1/3 V (2.20 b) 

00 

Eqn. (2.19) is called Betz efficiency and this does not represent 
the max. efficiency of the machine since the mass flow rate 

through the disc is not AV^ but AV. Hence, the efficiency i.e. 

power output divided by power available is given by -- 

P „ = 4a(l - a) .... (2.21) 

(1/2) p A V 

00 00 

The max. efficiency is 100% at a = 1/2 which yields a power 

coefficient (Gp) of 0.5. Thus the efficiency at the max. power 
coef f icient ( i.e. = 16/27) is 88.8% .( Putting a = 1/3 in 

Eqn. (2.21)). 

The flow pattern and the flow properties of wind along the 

horizontal axis of a wind turbine are shown in Fig. 2.2 . 

2.1.2 Definition Of Power and Thrust Coefficients : 

The conventional definition for thrust and power 
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coefficients are, as given before. 


;i,/2) p A V'^ 

^ 00 00 


Cr. = P 

(1/2) p A 

00 00 


Shankar [1976] has introduced two new coefficients Cp' and C,^' 
based on local wind speed at the rotor disc. He redefines the 
power coefficients C„ and C„' with a factor k = 16/27 in the 

denominator in order to normalize w.r.t. the theoretical max. But 
, we shall not use his definition but stick to convention. 
Therefore, we define 

C ^ = P , 

P (1/2) p A V 

His definition of and are according to convention, 

therefore 


Tl/2) p A 


. . . ( 2 . 22 ) 


Making use of Eqns.(2.17) and (2.18), 

(V^ /V)^ 

= 4 (V /V - 1) (using Eqn.(2.11)) 
00 

Hence, 


(V/V ) 
00 


1 + C^'/4 


(2.23) 



CHAPTER - 3 


3.1 Parameter identification : 

3.1.1 Rotor Area required : 

Any Wind Energy Conversion System (WECS) will have the follow! 
subsystems : 

Rotor, power transmission system and electric generator. 

Each of these subsystems has its own efficiency. Therefore, the 

total power extracted by a WECS may be written as 

P = Cp E^ E (1 / 2) p A (3.1) 

where , 

Cp = wind rotor power coefficient 
E = transmission system efficiency 

E^ = generator efficiency 
A = swept area required 
= wind speed 

= air density at the rotor hub. 

The existing wind turbines at Tuticorin are at about 20 m. 
elevation from the ground and hence the density at 20 m. altitude 
is not very different from sea - level density of air, p^ = 1.225 
kg/m . A modern and well designed wind turbine will have Cp 

around 0.45 . Also, we can choose the values of E and E as 

g ■ ■trans,,,,, 

0.85 and 0.75 respectively for practical purposes. From wind 
energy survey conducted in India, especially in Tuticorin we 
conclude that the optimum rated wind speed is approximately 10 m/s 
[20] . So the design rated wind speed could be taken as 10 m/s ^ ^ ^ ^ 
After choosing the above parameter.s, the rotor area required for 
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10 kW wind turbine can be found as 
P 


A = 


(1/ .2) p C^E E 

00 CO P trans g 


Hence, A “ 57 m 


57 m ( after calculation ) 
2 


3.1.2 Rotor radius. Blade length and chord : 

Sir Robert Mcalpine & Sons Ltd. [1986] designed a 25 m. , 130 
kW VGVAWT for Department of Energy, U.K. They arrived at some 
optimum ratio between the radius of the rotor and blade length 
after iterative calculations which are attributed to the 
structural consideration of the wind turbine. So it is reasonable 
to choose some proportions for the radius and blade length of the 
10 kW VGVAWT as they did. Accordingly, following values of the 

rotor radius and blade length have been chosen 

Rotor radius , R = 4 . 5 m 
Blade length = 6.5 m 

In the same manner the chord is also chosen as -- 
Chord, C = 0.45 m 

Hence, the rotor area A = 9 X 6.5 m^ ( refer Fig. 3.1) 

- 58.5 m^ 

3.2 General Arrangement of the VGVAWT : 

The general arrangement of the VGVAWT is as shown in Fig. 3.1 . 
It is proposed to have the no . of blades, N = 2 . The horizontal 
telescopic cross-arm supports the two half blades on one side. On 
the other side the cross-arm is connected via a spring to a point 
on the -rotor axis. Each half blade is Gonnected via a control 
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strut to the sleeve which houses the telescopic cross-arm {Fig 3.2) . 
The whole arrangement of rotor is supported on a cantilever tower 
(Fig. 3.1 ) . 

3.3 Airfoil Selection : 

The VAWT extracts power both from the wind upstream and 
downstream of the tower, hence the airfoil encounters not only 
positive angles of attack but negative angles of attack also. The 
performance of the rotor depends on the (L/D) of the airfoil since 
it uses the lift force to generate power. The cambered airfoil has 
poor (L/D) for negative angles of attack for a given conditions. 
Therefore we select a symmetrical airfoil which has a reasonably 
good (L/D) for both positive and negative angles of attack. 

NACA 0012 airfoil has been chosen for the computation. The 
blades are fixed at zero pitch to the circular path it traces. 
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4. 1 Formulation of the problem : 

The behaviour of a VGVAWT shown in Fig. (4.1) can be 
qualitatively understood by studying the figures(4.2) & (4.3). 

These figures show elemental section of a blade at various 
azimuthal position, (p, for a given wind speed V^and at a 
particular angular speed Q. When the local azimuthal speed rfi is 
large compared to the local wind speed V, the blade element is 
unstalled for all (p as the effective angle of attack a is small. 
In this condition, the lift dL contributes positively to the 
torque while elemental drag dD detracts from it. On the other 
hand when rQ < V, the effective angle of attack can vary from 0° 
to ± 180°, ie, the blade may be stalled and or in reversed flow 
over part of its trajectory. 

4.1.1 Expression for Radius and Area in terms of Arrowhead angle |3 

The expressions for R & A are found in terms of the reef 
angle T (see Figs. (3.2), and (4.4)) as follows; Since radius is 
changing with F, it is denoted as R' (Fig. 4.4) to differentiate 
from R which refers to the radius of the rotor (Fig. 3.1) 
corresponding to F = 0 and likewise swept area is denoted as A' , 
where A represents the swept area when the turbine is unreefed. 
As the strut links the blade (at C.G.) and the cross arm it is 
necessary to find its length for F = 0. In order to 
counterbalance as far as possible the centrifugal force acting on 
the rotor blade, the strut should be connected as near to axis of 
rotation as possible so that it makes a low angle. Since the tower 
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has diameter of the order of R/lO it is not practijilly possible 
to fix the strut exactly at the axis. Therefore, the strut is 
fixed at a distance R/lO from the axis of rotation. From the 
Fig. (4.4), the strut length, L for T = 0 (or = 90) can be 

calculated as 


L . /(^ R) = + 

where b is the quarter blade length. 


For b - 1.62 m and R = 4.5 m , L becomes 4.3 6 m and hence G works 
out to be 0.37. Refering to Fig. (4.4), one can write an 
expression for R' as 


R' = — yq + L COS0 + b cosp 


where 9 is called strut angle and /3 arrowhead angle. Note r=90°-/3. 
It is convenient to work in terms of the arrowhead angle. 

Also, L sin© = b sin/3 
i . e . 

cos© = }/ 1 - sin^/3 
Therefore , 


R' 


R 


10 


+ L 


G^sin^/3 + b cos^ 


From Fig . (4 . 6) , we have for /3 = 90°, 


0.9 R 


cos© 

o 


(4 .la) 


L 



b 

0. 9 R 


tan© and 
o 


(4.1b) 


— = G = Sin© (4.1c) 

. Li o 

where, © = ©^ when = 90° 

Hence, R' can be expressed as 

R' = -^TT + " 'il - sin % sin ^{3 + 0.9 R tan © cos (3 

10 cos © o o' 

o 

i.e. R' = ^ (4.2) 

where, f. = — ^ + 0.9 J sec^© - tan^© sin^/3 + 0.9 tan © cos (3 
1 10 o o o 

(4.3) 

The swept area A' can be found refering to Fig. (4.5) as 

A' = (2 R' X 4 b cos D - (8 b^cos F sin F) 

i.e. A' = 8R'b sin^ - 8b^ cosp sin^ 

i.e. A' = 8b sin^ (R' - bcos/3) 

Using Eqns.(4.1b) and (4.2), the above Eqn. becomes 

A' = 8 X 0.9R tan©^ sin^ (R f^- 0.9R tan©^cos/3 ) 

i.e. A' = 8R^ X 0.9 tan©^sin0 (f^ - 0.9 tan©^cos/3 ) (4.4) 

Also swept area corresponding to unreefed conditions is given by 


A = 8Rb^^^^ ^^ ^ ^ ^ ^ ^ 

Substituting for b from Eqn. (4.1b), we get 
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A = 8R^ X 0 . 9 tane (4 . 5) 

o 

After introducing Eqn. (4.5) into Eqn. (4.4), we get 
A' = A sin/3 [f^ - 0.9 tane^ cos/3] 
i.e. A' = A f 

where = sinp [f^ - 0.9 tan©^ cos/3] 

4.1.2 The Radius of a Blade Element: 

Let us consider an element ds at (r, - z) in the Fig. (4.4) . The 
analysis is first carried out for this lower half blade and then 
it can be generalised for the whole wind turbine. The radius of 
a blade element, can be written as 

r = R' + z cotfi 

where z varies from - 2b sin^ to 0 . 

By non - dimensionalising with the unreefed radius R, we get 

r = R' + z cot /3 (4.8) 

where r = ^ , R' 

K. 

4.1.3 Blade Element Analysis: 

We define the following unit vectors. 

e = unit vector along x-direction 

X 

By. = unit vector along y-direction 

e = unit vector along z-direction 

z 


R' 

R 


& z = 


R 


(4.6) 

(4 . 7 ) 
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e = unit vector along the chord of the airfoil blade section 
c 

e = unit vector along the normal to the blade axis and 

n 

perpendicular to the blade chord. 

The co-ordinate system is shown in Fig. (4.6) . It is obvious that 
e^ makes an angle (3 with the vertical. 

e^ = - cos^ e^ + sinp {cos<p e^ + s±n(p e^.) (4.9) 


e = - sind> e + cos(f> e,, 
c X y 


( 4 . 10 ) 


and for the elemental length ds located at a radius r from z-axis 


ds 


dz 

sin0 


also 


sin0 


1 



(dr/dz) 


2 


(4.11) 


Under unreefed conditions, the relative velocity at the blade 
section, is given by 

U = - V e + rn e - co3<f> e ) (4.12) 

R y ^ y 

However, the wind speed contributing to the lift and drag is given 
by the component of U^^ in the plane of the blade section. This 
component U is given by 
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ie, U = (rQ + V cos0)^ + (v sirup sin/3) 


(4.13) 


The effective angle of attack a at the blade section is given by 


tan a = 


^R-^n 

^R-^c 


V sin/3 sin (p 
rQ + V cosp 


(4.14) 


4.1,4 Equation of Mean Power and Mean Thrust: 

If the blade is of chord C with sectional lift and drag 
coefficients C^and C^, then the elemental lift and drag are given 
by 


dL = (1/2) pJJ^ C ds (4.15) 

dD = (1/2) pJJ^ C ds (4.16) 


At the azimuthal position the quasi-steay torque dQ and thrust 
dT due to element ds are given by 

dQ = r {dL sin a - dD cos a} 

dT = - cos <p {dL sin a - dD cos a} + sin /3 sin {dL cos a + dD 

sin a } 


respectively. 


By Eqns. 

dQ 


(4.15) and (4 .16) , 
r (|)p^U^C |c^ sin a - 


cosa}" ds 


} 
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dT = 

- coB(() (C^ sina - cosa) | ds 

Let sina - cosa and (4.16a) 

C = Ct cosa + Cj sina (4.16b) 

n 1 a 

where, and are tangential and normal force coefficients 
respectively. 

Hence , dQ = r(l/2)p C ds 

and, dT = (l/2)p C (sinp sin<^ - cos(p C^.) ds 

Using Eqn. (4.11) , the above eqns . become 

dQ . r(l/2)p u=“ c ^ 

and, dT = (-|- p C)(sin(3 sin0 

Therefore, the mean power and mean thrust produced by the lower 

half blade of Fig. (4.4) are found by integrating and averaging the 

above equations over interval 0 to 2ii and also integrating from 

z = -z = -2b sin/3 to z = 0. Thus 
1 



(hp ifc 

Z 00 


\ '' 

^ si 


inp sin<^ 


(Cj^cosa 


sina) 
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0 

p = - j- 

hb 271 

-z 

1 


2n 

tQ 

o 




d(j> dz 
sin/3 


and 


“hb 
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r)P, 


U^C 


sin/3 sin(p C^- 


cos<p 



d(p dz 
sin/3 


where the subscript hb signifies value for half -blade. Since each 
blade has two equal half blades, for the mean power and mean 
thrust of the whole wind turbine, the above expressions are to be 
multiplied by 2N, where N is the number of blades. Thus the mean 
power and mean thrust of a wind turbine can be expressed as 
follows : 


-z 

,.1 


2tt 


P = - 


N 


n 


rfl (hp U^C C. ^ - ,^1 and 
2 ^00 t sinfi 


(4.17) 


o o 


T = - 


N 

n 


- 2 271 

1 


O O 




sin0 sin(f> C^- cos<p 


d<^ dz 


sin/3 


(4 .18) 

If we define the solidity cr, the local tipspeed ratio A' , local 
power coeff.Cp' and local thrust coeff. as follows : 


cr 


N C 
R 


N C 
R' 


(4 . 19a) 
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f2 R Q R' 

V - 


(4 . 19b) 


C 


P 


P 

(1/2) p A'V^ 

^ 00 


(4 . 19c) 


C 


T 


~r 

(1/2) p A'V^ 
00 


(4 . 19d) 


It is to be noted that Cp' and C^' have been defined for reefed 
case based on the local speed V and therefore, area A' has been 
used for non - dimensionalization . Earlier, Selvaraj (1992) used 
the area A, i.e., unreefed intercepted area of wind turbine , for 
non-dimensionalization, i.e., he used the following definition of 
Cp' and C^' 


S' (l/2)p A 
00 


and C^' 


T 

(1/2) p A 
00 


But, as it is obvious that as the wind turbine reefs, the 
intercepted area decreases resulting in a decrease in absorbed 
power and hence decrease in power coefficient. Therefore, the area 
A' has been used for non-dimensionalization of Cp' and C,p' as in 
Eqns. (4.19c & 4 . 19d) to take into account this fact. Selvaraj 's 
analysis does not incorporate this fact because he had used area A 
for non-dimensionalization of Cp' and C,p' even in reefed case. 

Here C is the chord, and N is the number of blades . 
Normalising the rotor area by R, we have 



- K . (_^)2 

f „ ^ R' ’ 


(4.20) 
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CJsing the Egns . (4.19) & (4.20) in Eqn. (4.17) and dividing both 

sides of Eqn. (4.17) by (1/2 p A'V^), we get 

00 


P 

1/2 p A'V^ 
00 


-z 271 
1 

Nc r 

Tt 

o o 


r n 1/2 p 

00 t 

1/2 p A' 

00 


d0 dz 
sin/3 


or 


C ' 


ie, 


- z ^ 2ti 


NC 




r 

QR I 

U 

] ’ dp 

dz 

7 T R 




R 

V 1 

V 

V 

) sinfB 

R 


o o 


C 


» 

P 



2tt 


a-y I 

TT J 

o 


o 


r 


(U/V)^ 


d(t> dz 
sin /3 


(4.21) 


Similar deduction can be done for also and we finally arrive 

at the expression for C' as 


-Z 2n 
1 


C • 

T 


Trk 


cr(U/V)^ (sin(^ sin^ C 


n 


0 0 


vhere (U/V) ^ = (rA' + cos <p)^ + (sin (p sin ^) 


- cos<l> Cj. ) 


d# dz 
f^sTrip 


(4.22) 

(4.23) 
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Also Eqn. (4.14) , after using Eqn.(4.19) becomes, 


tan a 


sin p sin (p 
r X' + cos <p 


(4.24) 


The above expressions namely, Eqns.(4.21) & (4.22) will now be 

converted in terms of useful (true) tipspeed ratio {X) by using 
Eqn. (2.23). The detailed analysis was done by Shankar [1976], 
which avoids iterative procedure while calculating Cpand C,j,- 
Thus , 


X = X' (V/V ) = 


X' 

(1 + Cp'/4) 


(4.25) 


Cp = CpMA'/A) (V/VJ 


C ' f 
P 2 


(1 + C^'/4) 


(4 . 2< 


Cp = C,p'(A'/A) (V/VJ 


C,„' f., 

T 2 


(1 + C^V4) 


{4.2’ 


The mean torque coefficient, may be deduced from Cp as follows 


Q n ^ 

i/a p A^V 
'^00 00 


, Q 2 

1/2 p A^V 
•^00 00 


n R 


V R 
00 


ie Cp = Cq a. 


(4.28) 


Hence 


A 



where C, 


1/2 p A V R 
00 00 


{4 .29) 


We note that for blades symmetrically curved about an 
of radius R' with const, chord and of symmetrical section, 
Eqns . (4.21) & (4.22) reduce to 


-Z 


1 n 


C ' = 
P 


2aX' 
n k 


(U/V) ^ r C. - t 

t sin /3 f. 


(4.30) 


0 • 0 




T 


IV J J (U/V)= (sln^ sine - cos^ C^) 


0 0 


From Eqn. (4.1b) , we have 


- 5 - = 0.9 tan© 

R o 


(4 . 31 : 


Also, k = ^ 


= 7.2 tan© 


Therefore, Eqns. (4.30) and ( 4.31) become 


O' X‘ 


z n 
1 


3 . 671 tan© 


I I ^ IwV, 


0 0 


& C^' 


O' 


Z 7T 
1 


3. 6 n tan© 


r r (U/V)^'(sin(;> sinp C - cos(^ r 
J J n c. X .s i n 


0 0 


j!in0 


(4.33) 


where -z = — z|/R 
: 1 : 1 


SinS n ^ ^ „ 

p ' - = -1.8 tan© sxnp 

R o ' 


2 b 
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4.1.5 Aerodynamic Forces : 

Load distribution along the span of the blade is necessary for 
structural design of the blade. These forces are nothing but the 
tangential and normal forces acting on the blade. 

The tangential force per unit length of the blade is given by 

F (0) = 1/2 p C. C (4.34) 

C 00 t 

and the normal force per unit length of the blade is given by 

F (0) = 1/2 p \f C C (4.35) 

n 00 n 

Now Eqn. (4.34) can be written as 


F^(^) = 1/2 p„ C (4.36) 

Similarly, the normal force becomes 

Thus for a given freestream velocity, V^, if the parameter 

(U/V )^ is known then the aerodynamic forces can be evaluated. 

, 00 ■ 

The expression for the above parameter may be deduced as follows. 


= (u/v)' (v/vj' 


(4 . 38) 
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2 

where (U/V) is given by the Eqn.(4.23) and 




Note that the normal force is positive for os (p n and negative 
for (p ^ 2n . 
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Results and Discussion : 


The primary interest of the present work is to evaluate the 
aerodynamic characteristics of VGVAWT and the aerodynamic loads 
acting on the blade of VGVAWT. 

The VGVAWT has the following specifications : 


Elevation 

20 

m 

Rotor radius 

R 

= 4.5m 

Rotor area A 

= 

58 . 5 m^ 

Blade length 

= 

6.5 m 

Chord length 

(C) 

= 0.45 

No. of blades 

(N) = 2 

Solidity cr 

= 

0.2 

(V ) 

00 rated 

10 

m/ s 

(Cp) f H = 

F rated 

. H = 4 

rated 

0. 

42 


(Power) = 10 kW 

rated 

(fJ) = 85 rpm 

rated 

Airfoil used : NACA 0012 

The eqns . (4.32) and (4.33) have been solved numerically using a 

FORTRAN program which has been developed. The values of Cp ajid 
were obtained for various values of X' for a given value of /3 . 
Then the eqns. (4 . 25) , (4 . 26) , (4 . 27) and (4.28) are used to 
calculate X,Cp,C.j, and C^. NACA 0012 airfoil has been used for 

the computation. The section characteristics (C^ and C^) have been 
adopted from Ref . C5] and reproduced in Figs. (5.1a) and (5.1 b) . 
These characteristics were obtained experimentally for Re. No. = 
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0.3 X 10 . The Re. No. for the present case is 0.38 X 10*”, 

corresponding to = lo m/s and hence the use of the above data 
is justified. 

In Fig. 5. 2 , the variations of Cp,C^ and with A are 

presented for (3 = 90°. Cp is maximum at tipspeed ratio (A) = 4 and 

this value is 0.52 . For values of A less than 0.8, C„ is close to 

ir 

zero. After that it rises with A till A = 4.0 and then, it starts 
decreasing afterwards. At A = 0, has a non - zero value of 

0.125 and it increases with A. approaches to a value = 1.0 

asymptotically over the range of tipspeed ratios from 6 to 
infinity. Also, (C-,) = 0.15 occurs at A = 3.2, before 

Qmax Praax 

occurs . 

The Figs. 5. 3 ,5.4 ,5.5 and 5.6 show the variation of Cp, 

and (V/V ) with the arrowhead angles respectively. Also the 
00 

results of Selvaraj (1992) have been overlapped for comparison . It 

is to be noted that for lower tipspeed ratios (near A = 0) , the 

curves of C„,C_,C,., & V/V for the present case and for Selvaraj 
p 1 y 00 

overlap. For higher tipspeed ratios, Selvaraj 's analysis 
overpredicts the values of Cp,C,j,,CQ & compared to the present 

analysis. 

From Fig. 5. 3 , it is obvious that decrease in arrowhead angle, or 
in other words, increase in reef angle results in decrease in 
maximum power coefficient for both Selvaraj and the present case 
(Pandey, 1996) which is obvious because the decrease in arrowhead 
angle will decrease the intercepted area thereby reducing the 
absorbed power. The max. power coefficient occurs at a lower 
tipspeed ratio when the arrowhead angle is decreased for the 
present case. For a given arrowhead angle, the max. power 
coefficient occurs at a lower tipspeed ratio for the present case 
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than that for Selvaraj ' s case. The smaller the (3 the more 
flattened is the curve and Cp goes to zero at a higher tipspeed 
ratio than the rated tipspeed ratio (A) =4.0 for both the 

rated 

cases . 

The Fig. 5. 4 shows that has a value of 0.12 at A = 0 and it 
starts increasing with increasing tipspeed ratio for p = 90° for 
both the cases i. e. the present case and the Selvaraj ' s case. For 
other values of /3, also increases with increasing A .The value 
of at the zero tipspeed ratio decreases with decrease in 
arrowhead angle for both the cases 

Fig. 5. 5 shows that at A = 0, the starting torque coefficient Cq is 
0.026 for 13 = 90° and it starts decreasing with decrease in /3 for 
both the cases . 

Fig . 5 . 6 shows that the velocity ratio decreases with 
tipspeed ratio for a given (3 for both the cases. At zero tipspeed 
ratio, increases and approaches towards 1.0 as f3 is 
decreased. 

For |3 = 90° , (V/V ) is 0.97 and for /3 = 30° , it is 0.99 for the 

00 

present case as well as for Selvaraj.. For 13 = 90°, 
decreases much faster with A than for /3 = 3 0°. 

The effect of solidity on Cp is given in the Fig. 5. 7. The decrease 
in solidity shifts the curve towards right and also the curve 
becomes more and more flattened near the maximum. The increase in 
solidity shifts the (C„) towards lower tipspeed ratio There 
exists an optimum solidity ( between 0.2 and 0.3 ) for which the 
power coefficient is max. for /3= 90° . Shankar ' s calculations 1 non 
- uniform induced velocity approach ) have been superimposed for 
comparison. It is found that the present theory overpredicts the 


maximum power coefficient . 
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Fig. 5. 8 shows the angle of attack variation of the blade with 
azimuthal position for different tipspeed ratios { (S = 90° ) . It 

is obvious that for zero tipspeed ratio, the angle of attack, a, 
is same as (f) and that's why we get a straight line showing <p = a 
for A= 0. Note that the blades stall for tipspeed ratios from 0 to 
2.3, the stall angle being 14°. The variation of oc with varying 
from 180° to 360° is obtained by double reflection of this graph, 
first about vertical plane ( (^ = 180° ) and then about horizontal 
plane { a = 0° ) 

Figs . 5 . 9 and 5.10 show the variation of the tangential and normal 
force coefficients, C^. and C^, respectively as functions of 
azimuthal angle (</») for various tipspeed ratios ( for /3 = 90° ) . 
The results have been plotted for a single blade. For A = 0, the 
value of is -0.013 for the blade at 0 = 0 and 0.013 for the 

second blade at 0 = 180°. 

For computing the tangential force and normal force per unit 
length of the blade, we need to calculate the value of (U/V^) and 

hence, it is also computed and plotted in the Fig. 5. 11. Note that 

the aerodynamic forces acting on the blades are greatest for the 
unreefed condition, therefore Cj.,C^ and (U/V^) ^ are computed only 
for ^ = 90° . 

In 1983 , ^ Mathews experimentally obtained power coefficients of 25 
m VGVAWT for different reef angles . This has been compared with 
Selvaraj and the present analysis as shown in Fig. 5.12 . From 

Fig. 5. 12 , it is obvious that for |3 = 90°, the 25 m VGVAWT has 

= 0 . 45 at the tipspeed ratio 3 . 5 whereas the present and 
the SeLaraj's 20 m VGVAWT has = 0.52 at tipspeed ratio 4. 

This discrepancy may be due to the as.su.nption of uniform disc 
velocity (V) in the present theory. Hence, we can say that the 
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present theory overpredicts the mnx. power coo f P i c rent: . But, the 
results of present theory are more closer to Mathews' experimental 
results than that of Selvaraj's (Fig. 5.12) . Hence the present 
analysis gives more reliable results than that of Selvaraj ' s 
analysis (1992) . Note that for /3 = 90°, the present analysis and 
the Selvaraj's analysis give the same results. 

From Fig. 5. 7 , one can infer from the graphs obtained by Shankar 
[1976] that the uniform disc velocity theory (present case) 
overpredicts the max. power coefficient by about 10% when the 
solidity ratio is 0.2. 

Also we can observe that the increase in the reef angle (i. e. 
the decrease in the arrowhead angle) results in decrease in 
{Cp)^^x the cases (Mathews, Selvaraj and present case, 
Fig. 5. 12), but Cp goes to zero at a lower tip.speed ratio in the 
case of 25 m VGVAWT (Mathews) whereas Cp goes to zero at a higher 
tipspeed ratio, for the present case and for Selvaraj's case. This 
may be due to the fact that the 25 m VGVAWT reefs without 
extending radially outward , thereby producing more rapid decrease 
in swept, area which results in rapid decrease in Cp with X. But, 
the present case VGVAWT has a special feature of extending 
radially outward when it reefs. This feature decreases the swept 
area gradually and hence Cp becomes more and more flattened near 
the maximum as reef angle is increased. 
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Conclusions & Recommendations: 

1. The results of present theory (refer C_ vs X curve of Fig. 
5.12) are closer to Mathews' experimental results than that 
of Selvaraj ' s results. Hence, the present analysis improves 
over Selvaraj ' s analysis. 

2. Selvaraj 's analysis overpredicts the values of Cp, C^, and 
(V/V^) compared to the present analysis. 

3. The present analysis and the Selvaraj ' s analysis have assumed 
the blade to be having 2-D characteristics all along the 
span. The analysis can be repeated by considering the blade 
as a wing to improve the results . 

4. Present analysis and Selvaraj 's analysis have adopted the 
uniform induced velocity approach. Future work can be done 
by adopting non-uniform induced velocity approach to improve 


the results . 
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Summary of the work done on VGVAWT at I.I.T. Kanpur 

6.1 Introduction 

The potential of Vertical Axis Wind Turbine (VAWT) was recognised 
in the early seventies. Since then development of this wind 
turbine has witnessed many changes . The variable geometry VAWT was 
investigated and tried out in the UK. A 17m diameter 100 kW 
variable geometry VAWT was set up in the Isle of Scilly ( Mays et 
al 1987 ) . This turbine had an active control system. Active 
control systems depend upon transducers to sense wind speed, rpm, 
reef angle etc . whereas the passive control systems do their own 
sensing. These transducers must be very rugged, weatherproof and 
reliable and thus active control systems are expensive (also they 
tend to be quite sensitive to lightning and coastal regions in the 
tropical and sub - tropical parts of the globe are lightning - 
prone) . In late 1987, there was a failure in t)ie control system in 
the turbine of Scilly Isle which led to an overspeed situation 
causing damage ( Mays 1989, Lindley 1988/89 ) . 

The variable geometry rotor was given up in favour of the fixed 
geometry passive control (stall regulated) rotor ( Lindley, 1989 ) . 
In a VAWT the blade aerodynamic load alternates in direction 
(radially inward and outward) and hence the more severe loads in a 
stall regulated fixed geometry rotor would lead to greater 
fatigue. Ghosh (1990) predicted quite early that the fixed 
geometry rotor may not survive in a coastal area swept freqiiently 
by Gyclones and tornadoes . Wind speeds in a tornado range from 110 
to 550 km/hr and those in a cyclone from 180 to 280 km/hr ( Sachs, 
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1972 ) . As it so happens, the fixed geometry VAWT in the UK, a 
region which has milder storms, shed blades due to structural 
failure (Page, 1991 ) and the VAWT program there is now on hold. 

The variable geometry VAWT subjects its blades to 
lesser loads and enjoys certain advantages especially in the 
coastal regions of the tropical and sub-tropical parts of the 
world swept by cyclones. In these regions, a Wind Turbine ( WT ) 
harnesses a slightly inferior average wind compared to the 
temperate zone (near 40 latitude) but must survive the tremendous 
wind speeds of a cyclone (Hurricane, Typhoon etc.) . Ghosh (1990) 
proposed a new method of passive control system called Centrifugal 
Reef Control which is likely to be reliable and cost-effective for 
the variable geometry VAWT. A passive control system has the added 
advantage of not being susceptible to lightning strikes. 

A prototype design project for a variable geometry VAWT 
with centrifugal reef control has been underway at I.I.T. Kanpur 
since 1993. This project is assisted by a grant from the Ministry 
of Non- conventional Energy Sources (MNES) , Govt, of India. Various 
aspects of this design are presented here. 

6.2 Rotor Design 

6.2.1 Control and Aerodynamics 

The wind turbine rotor assembly, as proposed by Ghosh 
(1990), consists of two telescopic cross-arms, each of which 
supports two half blades on one end (Fig. 6.1, the half blades are 
attached to an end plate which in turn is hinged to the cross-arm 
head) . The other end of the cross-arm is connected to the rotor 
axis through a spring. Each half blade Is connected via a control 
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strut to the sleeve which houses the telescopic cross-arms . The 
reef control mechanism is based on the centrifugal action wherein 
the centrifugal force acting on the blade cross-arm assembly is 
conterbalanced by a horizontal spring restraint . Any increase in 
wind speed beyond the rated speed leads to an increase in the wind 
turbine rpm above the rated rpm. This will increase the 
centrifugal force on the blade cross-arm assembly. As a result, 
the spring restraint is overcome and the cross-arm moves radially 
outwards and the blades reef to resemble an arrowhead. The net 
effect is to reduce the intercepted area and the component of the 
wind in the plane of the blade section, and thereby control the 
power output . 

The preliminary aerodynamic investigation of a VAWT with 
centrifugal reef control to produce a power of 10 kW in 10 m/s 
wind has been carried out. Airscrew theory ( commonly applied for 
aircraft propellers ) in modified form has been used to compute 
power coefficient . Earlier, the power and thrust coefficients 
were computed for various reef positions (given by arrowhead angle 
0 in Fig. 6.1) by Selvaraj (1992) using unreefed intercepted area 
(corresponding to 0 = 90°) for non-dimensionalization . That means, 
he used power and thrust coefficients defined as follows - 

S' = (1/2) p A V"^ S' W2)p~Xv^ 

(Xt w 

where, P = power extracted from intercepted area 

T = thrust acting -on the wind turbine 

p = free stream air density 

' ' 00 , 

A = unreefed intercepted area of wind turbine 
(for (3 = 90°) . 

= air velocity at the wind turbine 


V 
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He used this definition even for f3 < 90 " and got the power 

coefficient variation with various arrowhead angles as shown in 
Fig. 5. 3. The analysis has been repeated by Pandey(1996) using 

reefed intercepted area (A' ) for non-dimensionalization of Cp' and 

C„' as follows - 
T 


C 


P 


P 

(l/2)p 

00 


and C^' 


T 

(1/2) p A' 

CO 


where, A' = reefed intercepted area of wind turbine 
(for < 90°) . 

Obviously, for {3 = 90°, A' = A. 

The above definition results in greater values of for any 

given 8 and thus, by virtue of Eqn.(2.23), lesser values of V/V 

00 

This in turn leads to smaller power. 

Using this definition of power coefficient, the variation of power 
coefficient with different arrowhead angles has been evaluated, 
using a Fortran Code, and the results are shown in Fig. 5. 3 . On 

comparison , we see that the analysis of Selvaraj slightly 
overestimates the power coefficient. 


6.2.2 Structure and mechanical configuration 

The design of such a wind turbine requires an appropriate 
mass distribution such that the telescopic cross-arm is much 
heavier (about 2.5 times) than the rotor half blade (Ghosh ,1990} . 

It further involves estimation of centrifugal force and an 

assesment of the spring stiffness capable of withstanding the 
centrifugal force for different reef angles. 

The major dimensions of the wind turbine have been shown in 
Fig. 3.1 . A detailed design study is carried out to determine the 
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size and relative position of various components of the rotor 
assembly. Fig. 6. 2 shows various features such as the bearing pads, 
spring, sleeve, cross-arm, end plate, blades etc. 

One end of the control strut is hinged at the C G of the 
half blade. The sleeve, which houses the cross-arm and spring, has 
a box beam structure with taper over the inner radius portion. The 
cross-arm weighing approximately 55 kg is made of commercially 
available standard aluminium tubes. It is given a tubular 
structure to facilitate support by boundary lubricated bearing 
pads made of cast iron. Since aluminium offers a poor bearing 
surface, it is to be coated with tin through a chemical process 
(Horn, 1967) . Fig. 6. 3 shows the top view of the cross-arm which 
consists of two aluminium tubes placed side by side. Also the 
figure shows the Z - link which is connected onto the cross-arm by 
means of a fork [inset of Fig. 6. 3 ]. The purpose of the Z - link 
is to limit to and fro motion caused by the reversal of 
aerodynamic load on each blade and help maintain symmetry about 
the axis of rotation. The rotor half blade (NACA 0012 airfoil 
section) which is 3.25 m long, has a chord of 0.45 m and a 
rectangular planform. It weighs approx. 22 kg. The .blade has an 
FRP composite structure, consisting of a skin wrapped over a D 
spar, meant to withstand fatigue (Fig.6.4) . It is provided with 
discretized internal taper to reduce weight and bring C G closer 
to the end plate (Fig.6.4) . It carries a max. normal and 
tangential aerodynamic loads of 480 N/m and 80 N/m respectively. 
The detailed load calculations have been carried out by Selvaraj 
(1992) for the unreefed case of operation of wind turbine because 
the aerodynamic loads would be maximum for unreefed case and 
therefore, other reefed conditions would be automatically taken 
care of. Fig. 6. 5 shows the variation of the centrifugal force per 
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unit mass of the half blade with cross-arm movement as it reefs. 
The location of C G of cross-arm (denoted by r in Fig. 6. 6 ) for 
the unreefed case is 1.621 m from the axis of rotation (Ref. 33) . 

Fig. 6. 5 shows clearly that the centrifugal force is only 
slightly non-linear with r. Analysis reveals that for the present 
case, a stiffness of 0.4 kgf/mm is required for each spring in the 
pair to withstand the centrifugal force for different reef angles. 
The max. displacement of the cross-arm between unreefed and fully 
reefed conditions is 1.7 m. The springs are stretched by 0.8 m, 
even in the unreefed position. Hence, they have to undergo a total 
elongation of 2.5 m. Such special springs are made by a firm 
called Stumpp Schuele & Somappa Ltd. of Bangalore. The chosen 
spring has a free length of 3.0 m and weight of 58 kg. The 
elongation in the spring due to the centrifugal forces generated 
by its own mass is approximately 8.0 mm. However, this elongation 
may increase from 8.0 mm to about 16.0 mm in fully reefed 
condition (see Appendix 1, taken from ref. 32) . 

Now, we can summarise the wind turbine specifications as : 
Unreefed radius R = 4.5 m 
Rated tip speed ratio A = 4 

Rated angular velocity = 85rpin 
Half blade length L^= 3.25 m 
Blade chord length = 0.45 ra 
Rated wind speed = 10 m/s 
Cut in wind speed = 4 m/s 

The wind turbine has no cut out wind speed, since at extremely 
high wind conditions the blades can reef by almost 90 deg., and 
still produce power. Other wind turbine characteristics are : 

Mass of the half blade M^= 22 kg 
C.G. location of the half blade (Fig. 6. 6) cg^^ 


1 . 3 ru 
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Width of the end plate W = 150 mm 
Thickness of the end plate t = 15 mm 

ep 

Length of the end plate L = 0.45 m 

ep 

Mass of the end plate M = 1.0 ka 

ep * y 

width of the sleeve W = 175 mm 

S 1 

Thickness of the sleeve t = 10 mm 

S 1 

Height of bearing sleeve or pad h = 20 mm 

S 1 

Space between the bearings , bear_s = 1.5 m 
Location of the strut hinge X = 2.0 m 

h 

Location of the near end of X-arm (Unreefed) X =150 mm 

X 

Mass of the X - arm M = 55 kq 

The following work was done by J. Singh (Ref. 33) and the same is 
being presented here. 

6.3 Dynamics & Response 

The purpose of dynamic analysis is to study the transient 
response of the reefing blades of wind turbine and see how quickly 
they settle down to their new equilibrium position at increased 
wind speeds. The approach is to determine the changes taking place 
along the blade's trajectory. To evaluate the increased rpm and 
formulate the dynamic equations, we trace the motion of the wind 
turbine blade, starting at some value of the azimuth angle ^ . 

Let the wind turbine blade be at station 0 at time t = 0 when 

the azimuth angle ^ (Fig. 6 .8) . The subscript 0 is used to 

denote quantities at station 0 at t = 0 . The wind velocity at 

the actuator disc can be obtained from the ratio (V/V^ ) which is 

a function of arrowhead angle p and the tipspeed ratio X (refer to 
Fig. 5. 6) . Let us assume that there is a sudden jump in the free 
stream wind velocity at station 0 from 10 m/s to 11 m/s (i.e. unit 
step jump) . This increase in free stream wind velocity will result 
in a corresponding increase in the wind velocity at the disc. 
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Knowing (3 & \ , this increased V at station 0 ( denoted by ) at 
the disc can be evaluated in the following way {from Fig. 5. 6) : 

(These values have been taken from Selvaraj ' s work (1992) for 
reefed as well as unreefed case of wind turbine operation because 
his analysis would give conservative structural as well as 
d.yn3.rnxc3.1 d-BSigri a.nd h6ncG t.h.0SG valuGs havG bGGn takGn) 

For = 10 ra/s, X = w R /V = 4.0, /3 = 90° gives /V = 0.72 

” ( 1 ) 

For V = 11 m/s, X = w R /V = 3.64,/3 = 90° gives Vj /V = 0.74 

00 O 00 00 

( 2 ) 

where, subscript j denotes the variable after the jump. The ratio 

Vj/V = 0.74 is assumed to be constant for the calculations 

00 

performed at all the subsequent values of azimuth angles and 

arrowhead angles. The increase from Vo to Vj will modify the 

velocity triangle at station 0 and consequently, the angle of 

attack increases from ao to aj (enlarged view of inset in Fig. 

6.8 ) and there is a corresponding increase in the aerodynamic 

torque from Q to Q . At rated wind speed conditions, the average 
O J 

aerodynamic torque over a cycle is equal to the generator torque 

Q and the blades rotate at a constant rpm. With the jump in the 

aerodynamic torque from to , the net increase given by 

“ Q ) causes the wind turbine to accelerate to increased angular 
g 

velocity. 

To determine the aerodynamic torque produced by the half 

blade, we divide the half blade into n number of equal segments, 
with the i*'*' segment being located at a distance 1^ from the end 
plate and r^ from the axis of rotation (Fig.6. 9) . 

Then, 
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Q.b ■ ,5, r_(c^)^ ,3) 

where, C^ = tangential force coefficient 
Ppo = free stream density 
= half blade length 
C = blade chord 

D 

Uj = relative velocity of the i*^^ blade segment in the 
plane of the blade section. 

The aerodynamic torque Q produced by a single wind turbine blade 
is given by 

Q = (4) 


Hence, Q - p_^(L^/N)C^ , 

where, r^ = R' - l^Cos^ (6 


reefed wind turbine radius (Fig. 6.7; 
+ L Cos© + eg Cos0 ... (7) 

h s s * 


where, © = Sin'^ [{w Cos/3/2 + (eg + t ) Sin/3 + h }/L ] 

s ep b ep z'^s 

The distance 1^ of the i^^ segment from the end plate is 


1 = (L /N) (i - 1/2) ..... (9) 

ID 

The wind speed contributing to the torque Q is given by 

Uj = (r^w + V Cosip)^ + (V Sin0 Sin/3) .....(10) 

The tangential force coefficient (C^) for the i’^*' blade segment is 
given by -- 

(C ). = c Sina - C Cosa. . . .(11) 

t. 1 11 d I 


. . . (11] 


where, a = tan~ [V Sln/3 Sin4> / (r^w + V Cos<^) ] 


C and C are the section lift and drag coefficients. 

1 d 

The aerodynamic torque and for the first blade at station 

0 can now be calculated by the use of equations ( 4 - 12 ) as 


Q ■ = p (L^/N) C^ S [U^] r [(C ) ] 

o 1 '^oa b b 1 o 1 t 1 » 

, : 1=1 
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Similarly, the aerodynamln torquea and Q^^for the second blade 
can be obtained by computing „_,u, and (c^)^ for 180 + ^ . The 

total aerodynamic torques and produced by the wind turbine 
are then given by -- 

Q = Q + Q „ 

o ol o2 

Q = Q + Q 

Jl J2 

The net increase in the aerodynamic torque at station 0 which 
causes the wind turbine to accelerate to an increased rpm is given 
by --- 

AQ = Q - Q (15) 

J 9 

If I denotes the sum total of the moment of inertia of all the 
components of the wind turbine rotor assembly, the angular 
acceleration 6 is given by the following relation -- 


H 

10 + 2 [ (2V w) Mr'] +2 [4w Z (V - 1 f CosD m r ] = AQ ... (16) 

XX 1 = 1*^ ^ ^ 

The last two terms in the above expression arise due to the 
Coriolis effect of rotation. 

Here,. = radial velocity of the X - arm 

w = angular velocity 

Ij = distance measured from the end plate (Fig. 6. 10) 

The Coriolis component of the acceleration for the X - arm is 
2V w. If M is the mass of the X - arm and r' the location of the 

X X 

c.g. of the X - arm (r' = r + T, where r is the c.g. location in 
unreefed condition and T is the travel or displacement of the X - 
arm; Pig. 6 .7) , the expression within square brackets in the second 
term of above eqn. represents the negative torque arising out of 
the Coriolis acceleration of the X^- arm, and hence it is included 
in the LHS of the eqn. This expression has been multiplied by a 
factor of 2 to account for. both the X-arms . 
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The expression within square brackets in the third terra of 
above eqn. represents the negative torque due to the Coriolis 
effect obtained for a single wind turbine blade. This term is more 
complex since the blade hinged to the end plate is moving outward 

with velocity and at the same time it reefs through an angle T. 

« 

If r is the angular velocity with which the blade is reefing^ then 

♦ 

~ CosF) gives the net velocity of the i^*’ segment in the 

radially outward direction. If r^ represent the mass and the 

distance from the axis of rotation of the segment 

respectively, then the torque for a blade arising due to the 
Coriolis effect of rotation is given by 

N 

40) Z (V -if CosD m r 

1 = 1 X 1 I i 

This term has been multiplied by a factor of 2 in eqn. (16) to 
take into account both the blades of the Wind Turbine. 

The moment of inertia I is the summation of the moment of 
inertia of all the rotating components of the wind turbine and is 
given by -- 

I = I +I +1, +I . +1^,^ .... (17) 

blade X-arm sleeve spring end-plate 

For the angular velocity, o) = w^,/3 = 90°,r = 0,f = 0 and V^= 0, 

eqns. (16) and (17) can be used to determine the angular 

acceleration 0 at station 0 at time t = 0 . 

Let us assume that the blade reaches station 1 after time At = 
0.1 sec (Fig. 6. 8). Then, the angular velocity at station 1 can be 
obtained as -- 
u = 0 At + w .... (18) 

1 O 

The new azimuth angle at station 1 is given by 
= { (Wj+ td^)/2}At + (/>^ . . . . . (19) 

The angular velocity increases from at station 0 to w^at 

station 1 . This causes an increase in the centrifugal force 

- arm in motion. Obviously, the reefing 


thereby setting the X 
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comes into play only from station 1 onwards . Hence , at station 1, 
the arrowhead angle (3 is still equal to 90°. For 13 = 90° and = 

(p^, eqns. (4 - 12) can be used to evaluate the torque at 

station 1. 

Now, let us assume that the blade takes time At to move from 
station 1 to station 2 (Fig. 6. 11). The increased aerodynamic 
torque which causes the angular acceleration from station 1 
onwards is given by -- 
AQ = Q, - Q .... (20) 

1 <3 

where, = aerodynamic torque at stn. 1 
= generator torque 

For the value of AQ given by eqn. (20) , (3 = 90 deg. and V = 0, 

X 

the angular acceleration 0 at station 1 can be calculated by 
eqn. (16) . Since the wind turbine blades are reefing, there is a 
continuous change in the aerodynamic torque and angular 
acceleration from stn.l onwards. To evaluate this variation in 
torque and angular acceleration during the blade's motion from 
stn.l to stn.2, we divide At (Fig. 6. 11) into small time intervals 
each of magnitude 5t (assume 6t = 0 . Olsec). Obviously, the blade 
starting from stn.l reaches stn.l_l in time 5t. In this time St, 
the blades reef through a certain angle T and the c.g. of the X - 
arm experiences a displacement x in the radial outward direction. 
We have to obtain the values of F and x at station 1__1 which will 
subsequently yield the value of the torque at station 1_1. 

The change in the aerodynamic torque AQ (AQ = Q^ ^ - Q^) can then 
be evaluated and used to calculate the value of the angular 
acceleration 9 at station 1 _ 1 . Starting from station 1, the 
Various steps involved in arriving at the values of AQ^^^and O at 
station l 1 are similar as involved in going from stn.O to stn.l. 
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we keep on repeating this process from stn. to stn. till the 
arrow head angle 0 stabilises to some value for all the following 
stations and the difference in the rated and the final angular 

velocity (corresponding to the new equilibrium value of S 1 is not 
more than 1%. 


Based on the above formulation, the response of the wind 
turbine rotor assembly to a step jump in the wind velocity was 
observed. We took Selvaraj's result for load calculation even for 0 
<90 because it would give conservative dynamical prediction. The 
free stream velocity was assumed to increase from 10 m/s to 11 m/s 
and the time history of m and /3 was plotted. Fig. 6. 13 (a) shows w 
(in rad/s) plotted against time (in sec.) . The cyclic variation in 
u) observed in the graph is due to the variation in torque along 
the path traced by the turbine blades . It is observed that the 
time period of the cyclic variation in w decreases near the peak 
to increased angular velocity. Fig. 6. 13(b) shows the value of 
rpm averaged over one rotation plotted against time. The mean rpm 
is observed to increase and then smoothly decrease to the rated 
rpm in a short span of time. Fig. 6. 14 shows the time history of 
arrowhead angle 0. The response is observed to be critically 
damped and the arrowhead angle settles down to its new equilibrium 
value with practically no oscillations. 

6.4 Grid connection and Starting 

The turbine turns an induction generator which is connected to the 
grid. The principal advantage of the induction generator is its 
ability to draw power from tlie grlcr t r j .ni nr t: l.ho Lurl>ino , which 
has a poor starting torque like most other VAWTs. When a wind 
speed of 4 m/s (equal to the cut in wind speed) persists for more 
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than 10 sec., a wind speed sensor activates the starting 
procedure. As told earlier, the WT has no cut out wind speed since 
at extremely high wind conditions the blades can reef by almost 
90 degrees and still produce power. 
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Fig.3-2 Details of Rotor 
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5-12 Comparison of power coefficient with experiments 
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Fig. (6.1) Reef Control Mechanism 



Fig. (6 .2) 


of the Rotor ; Front View 
(Z - Link omitted for clarity) 
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Fig. (6.14) Time histoiry of arrowhead angle {(3) 
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APPENDIX 1 : Elongation of Rotating Spring 

Total number of coils •= 220 
Free length “ 3000 mm 

Therefore, eacli coil length = 3000/ 220 == 13.6 mm 

Mass of the spring = 58 Kg 
Mass M of each coil = 58/220 = 0.26 Kg 
Spring stifness = 0.4 Kgf/min 

Stiffness of each coil = 0.4x9.81x10^x220 == 863.28x10" N/m 

Centrifugal force in the coil — Mrw^ = 0.26 xrxw^ N 

where r= 13,6xl0~" m and u) = 8.9 rad/s 

Elongation Ar of the coil = 9.i262Lr2L!ilf 

863.28x10^ 

Ar == 3.01 176 X 10 rw^ meters 
Centrifugal force due to the 2"*^ coil = M(2r)w^ Newtons 

Therefore, elongation of the 2"^ coil = 2Ar and net elongation upto the 2"^^ coil = 
2 A r *4- A r 

Similarly, net elongation upto the 3*^^^ coil = 3Ar+2Ar-|-Ar 

Therefore, net elongation upto the 220^^ coil = 220 A r + ... + 3Ar+2Ar4-Ar 

Hence, net elongation in the spring due to the centrifugal force generated by its own mass = 

|220^21j^j. _ 24310 Ar 

= 7.88 X 10 meters 
8.0 mm 

When fully reefed, the spring stretches approximately to twice its length. I'he value of r in 
that case will approximately be equal to 27 mm and the elongation in the spring (rpm remains 
constant) due to its own mass goes upto 16 mm. 
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